Role of interactions in ferrofluid thermal ratchets 
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\^ • Orientational fluctuations of colloidal particles with magnetic moments may be rectifled with the 

' help of external magnetic fields with suitably chosen time dependence. As a result a noise-driven 

, rotation of particles occurs giving rise to a macroscopic torque per volume of the carrier liquid. 

. We investigate the infiuence of mutual interactions between the particles on this ratchet effect by 

_ ' studying a model system with mean-field interactions. The stochastic dynamics may be described 

by a nonlinear Fokker-Planck equation for the collective orientation of the particles which we solve 
approximately by using the effective field method. We determine an interval for the ratio between 
coupling strength and noise intensity for which a self-sustained rectification of fiuctuations becomes 
possible. The ratchet effect then operates under conditions for which it were impossible in the 
, absence of interactions. 
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PACS numbers: 82.70.-y, 75.50Mm, 5.40.-a, 5.60.-k 



INTRODUCTION 



, Thermal ratchets or Brownian motors are devices that are able to rectify thermal fluctuations and thereby extract 
directed motion from irregular microscopic chaos. Banished from equilibrium by the second law of thermodynamics 
they may nevertheless be found in various non-equilibrium situations (for a review see 01 )■ Having being used for 
a long time as thought models to highlight some subtle points in the foundations of statistical mechanics j3| they 
' became recently objects of more practical interest due to their possible relevance for biological transport mechanisms 
, 0,01 and potential applications in nano-technology 
Q ■ Among the many devices that have been proposed as Brownian motors j3| those based on ferrofluids 0, 0| are of 
O , particular interest since they allow a rather direct observation of the ratchet effect on the macroscopic level. Ferroflu- 
ids are suspensions of ferromagnetic nano-particles in carrier liquids like water or oil combining the hydrodynamic 
[ behaviour of Newtonian fluids with the magnetic properties of superparamagnets ^3|. The orientational Brownian 
^ ■ motion of the ferromagnetic grains can be rectified with the help of a suitably chosen time-dependent magnetic field. 
' Due to the viscous coupling between the rotation of the magnetic grains and the local vorticity of the hydrodynamic 
, flow the angular momentum produced by the ratchet is transferred to the carrier liquid and may be detected as 
hydrodynamic torque per fluid volume 

Although the interactions between the ferromagnetic particles mediated by the carrier liquid are of crucial impor- 
' tance for the macroscopic manifestation of the ratchet effect it was neglected in the theoretical modeling done so far 
! 1^- III the present paper we want to remedy of this drawback and to elucidate the influence of direct (dipole- 
■ dipole) and indirect (hydrodynamical) interactions between the ferrofluid particles on the ratchet effect in ferrofluids. 
a [ This will be done using a simplifled mean-fleld model of interactions which allows to analytically discuss some of the 
^ . qualitative changes that occur. 

I ' General aspects of the interplay between the ratchet effect and interactions were already discussed in [13 and P|. 

, In particular it was shown that under certain conditions the interactions may give rise to a spontaneous breakdown 
g ■ of symmetries that formerly inhibited the ratchet effect. The relevance of interactions in biological applications was 
Q ] investigated in 01 • 

Here we will flnd by studying a concrete, experimentally accessible example of a ratchet that depending on the ratio 
between coupling strength and noise intensity qualitatively different regimes are possible. In some of these regimes a 
self-sustained ratchet effect may occur under circumstances for which no such effect would show up in the absence of 
' interactions. The discussion will be centered around a non-linear Fokker-Planck equation for the stochastic dynamics 
5^ , of the collective orientation of the ferrofluid particles J3 which is approximately solved by adapting the well-known 
effective-fleld method from the theory of ferrofluids to the analytical investigation of the ratchet effect in these 
systems. 
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The paper is organized as follows. In section II we introduce the model and the basic equations and derive the 
non-linear Fokker-Planck equation for the collective orientation of the ferromagnetic particles. Section III contains 
some numerical solutions of this Fokker-Planck equation providing first insight into the different regimes that are 
possible in the system. In section IV we introduce the effective-field method and show how it may be adapted to 
the present situation. Using this method we will characterize in Section V the different regimes found in section III, 
determine their stability and discuss the different manifestations of the ratchet effect. Finally, section VI contains our 
conclusions. 



II. THE MODEL 



We consider the overdamped dynamics of N identical spherical particles with magnetic moments m.i,i ~ 1,...,N, 
dispersed in a Newtonian liquid with viscosity 77. The particles are under the infiuence of a time-dependent external 
magnetic field with constant a:-component and time-periodic y-component 

U={H,,Hy{t),0) , Hy(^t+^^ =Hyit) . (l) 

A suitable example for the time dependence of the y-component is 

Hy(t) = acos{nt) + Pcos{2nt + 5) , (2) 

with the scalar parameters a, /?, and 6. 

The change of orientation Ui = m.i/mi of particle i is given by 

dtUi = uJiX Ui , (3) 

where denotes the angular velocity the particle. 

The particle orientations tend to ahgn with the magnetic field. Since the latter is confined to the x-y plane the 
average orientations will have zero z-components and the main effects will show up in the time-dependence of their 
X- and y-components. To keep the analysis simple we will therefore only consider the dynamics in the x-y plane and 
put Ui^z = for all i. It is then convenient to parametrize the orientations in the form 

(cos 
sin^, . (4) 
J 

In the overdamped limit the angular velocities uji are determined by the balance of torques 

7nagn,i ^" ■N.j^,2'sc,?' ^ stoch,i ^" "^int^i • (5) 

The different contributions in this equation denote the magnetic torque due to the interaction with the external 
magnetic field 

N,„agn,j = TO Mi X H = m[Hx sin 0j - Hy{t) COS 4>i)ez , (6) 
the viscous and stochastic torques 

N„sc/t = -er/Vwi = -6T]Vdt(t) Bz , (7) 

and 

^^stoch,^ = yi2^V^6(0 ez , (8) 

respectively describing the interaction with the carrier fiuid, and the torque arising from the interaction between the 
particles 

K ^ K ^ 

Ni„t,i = ]^ X! ^ TV ^ sin(0,; - (j)j)ez ■ (9) 
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In the above expressions denotes Boltzmann's constant, T temperature, and V the volume of the particles. The 
^i(t) are identical and independent Gaussian noise sources with zero mean and correlation 

Um,it'))=5^,5{t~t'). (10) 

The interaction term (jsj is of mean-field type as e.g. also advocated in Although the direct magnetic dipole- 

dipole interaction as well as the indirect hydrodynamic interaction between the particles are much more compHcated 
than the simple assumption (j^l several interesting implications of interactions will become apparent already in this 
mean-field description. 

>From Q we find the following set of Langevin equations describing the stochastic dynamics of the system 



^ [mHj, + j^^cos <j)j^ sin <j)i - (mHy{t) + ^-j) + V2D £.i{t) (11) 



6riV\ N j-^^ "V 6i]V 



with the diffusion coefficient D given by 



D:^^. (12) 



It is convenient to introduce dimensionless quantities by performing the rescalings t t/fl, H G'qVil./m'H., 
D n'^D and K 6i]Vn K to obtain 



dt(l>i = (^Hj; + ■^Ycos(l)j'^ sin0j - (^Hy{t) + -^X^^^^'^j) '^°^4>i + V2D ^iit) . (13) 

Equivalent to this set of Langevin equations is a Fokker-Planck equation for the joint probability distribution W{(l>i,t) 
of the particle orientations. 

As usual in a mean-field model we are interested in the limit N ^ oo. It is then useful to introduce the distribution 
function for particle orientations 



1 " 



i=l 



As discussed in detail in [l^ IT^ IT^ P{(1>, t) becomes for N oo independent of the specific realization of the 
noise d (t) (self-averaging property) . Consequently the same holds true for the collective orientation of the particles 
defined by 



Sit) ^ J^J2Mt) = dc^uP{c^,t) = I 



cos ( 

= I (sin0) I , (15) 

where (.) denotes the average with P{(j),t). We therefore get from Ijl3|l decoupled Langevin equations of the form 

dtA = {H, + KS,{t)) sin0, - {Hy{t) + KSy{t)) cos0, + V2l5 ^.,{t) . (16) 

The solution W{4>i,t) of the equivalent Fokker-Planck equation now factorizes, W{4>i,t) = HiLi w{(j)i,t), with w{(j)i,t) 
obeying the single-particle Fokker-Planck equation [l^ [T^l 

dtw{^, t) = (u.(0, t)d^U{dp, t)^ + D dlw{<p, t) . (17) 

The potential U is defined by 

C/(0, t) = -u- (H(t) + KS{t)) = -{H, + KS4t)) cos (f> - {Hy{t) + KSy{t)) sin0 , (18) 

and depends parametrically on H(i) and S{t). 

We may now average ifTljl over the realizations of the (t)i [t) to find that P obeys the same equation as w, 

dtPicj), t) = a4P(0, t) d^UicI,, t; U{t), S{t))) + D dlP{^, t) . (19) 
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FIG. 1: Modulus S of the collective orientation S defined in l|15|l as function of the ratio K/D between interaction and disorder 
strength. In the disordered phase, K/D < 2, 5 is identically zero, for K/D > 2 a spontaneous ordering takes place in the 
mean-field model. 



This equation is closed by ltT5|l . As characteristic for a mean-field system we have thus reduced the dynamics of 
N ^ CO coupled degrees of freedom to the dynamics of a single degree of freedom in an augmented field 



to be determined self- consistently. For a related application of mean-field techniques to ferrofiuids see also j2£ 
For a static field, H(t) = H, eq. Ijl9|l admits the stationary solution 



1 



exp 



mf^x 



C0S( 



mf,y 



sm( 



D 



(20) 



(21) 



D 



Here and in the following /„(a;) denotes the n-th Bessel function of complex arguments. In the absence of an external 
field, H = 0, the self-consistency condition ((T5ll takes the form 



(sin< 



= B 



KS 

IT 



s 



where we introduced the function 



B{x) 



h{x) 



which is a monotonously increasing and satisfies 

B{x) 



X 

2 

■ 1 



for 
for 



oo 



(22) 
(23) 

(24) 
(25) 



Eq. (l22|l coincides with the self-consistent equation for the ferromagnetic mean-field x-y model. The dependence of 
the modulus S of the collective orientation S on the ratio between interaction and fiuctuation strengths is shown in 
Fig. ^ Note that the direction of S in the ordered phase is arbitrary. Note also that for short range interactions the 
situation is rather different (21I . 

The central quantity of interest in connection with the ratchet effect in the present system is the time and ensemble 
averaged torque transferred from the magnetic field to the particles 0, 1^ 



(N) = lim - 



dt{u) x'H. = - lim — dt {dtcp) , 



(26) 
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where the second equality follows from Ijlfill . 

I n the absence of interactions between the particles the necessary conditions for the ratchet effect to operate, i.e. 
for (Nz) ^ 0, have been discussed in detail in [aE^- In particular it was shown by a symmetry argument that = 
implies (N^) = 0. In the present mean-field model the role of H is played by H„y = H + S' and it might therefore 
be possible to find (Nz) ^ whenever Hmf,x 7^ 0. However, the latter condition may be fulfilled by = and 
Sx ^ 0. Hence we may suspect that rectification of fiuctuations may take place in the interacting system even under 
conditions for which it would be impossible in a system without interactions. Note that there is no trivial mapping 
between the two cases since Umf.x is time dependent and has to be determined self-consistently. Nevertheless the 
detailed investigations discussed below will show that the interactions between the particles may indeed give rise to 
a non-zero value of Sx which in turn may drive the ratchet effect in the system even if Hx = 0. 



III. NUMERICAL SOLUTION OF THE FOKKER-PLANCK EQUATION 

In the present section we investigate some general features of the dynamics of the system under consideration with 
the help of a numerical solution of our central equation (I19II . To this end we expand P{4>,t) in Fourier modes with 
respect to (j) 

00 

exp(m(/)) (27) 

n— — 00 

with time-dependent complex expansion coefficients a„(t). From (tTfijl we then have 

Sx{t) ^2T:^ai{t) (28) 
Sy{t) = -2T:'iai{t) , (29) 

whereas the average torque is given by 



{Nz) = iJx^ai + Hy^ai (30) 

Using the Fourier expansion of -P(0, t) in (I19|l we obtain an infinite system of coupled ordinary differential equations 
for the coefficients a„ [t) of the form 



where 



Starting with 



{dt + Dn^)an{t) = -{g{t)an-i - 5*(t)a„+i) . (31) 



g{t)=H^,-iHy + 2TTKai{t). (32) 



0-0 = TT (33) 



as required by normalization of P{4>, t) we may solve the hierarchy of equations Ij31|l iteratively up to some value Umax 
of n. Using different values for Umax we have found that rimax = 10 is sufficient to get accurate numerical results. 

As discussed at the end of the previous section the case Hx = Q is of particular interest. Depending on the ratio 
between coupling strength and noise intensity we find in this case three different regimes which are characterized by 
the time dependence of the average orientation S{t) as shown in Fig.|2| In the first regime, K/D < 2, we find Sx{t) = 
and Sy{t) changes sign during one period of Hy[t). In the second regime, 2 < K/D < 3.7, Sx acquires non-zero values 
whereas the behaviour of Sy{t) is qualitatively similar to the first case. In the third regime, 3.7 < K/D, we have 
again Sx = but now Sy{t) is pos itive for all t. Similar regimes have been discussed also for the spherical model j23| 
and the anisotropic x-y model [2^ in time dependent external fields. 

As discussed in the previous section the symmetry analysis performed in Q implies that for ^ and Sx = no 
ratchet effect is possible. Accordingly, using (|3(11l we find (Nz) = in regimes 1 and 3. In regime 2 we have Sx 
which allows non-zero values of {Nz). This is indeed what we find with, however, two peculiarities. First, the values 
for (Nz) are extremely small as long as K/D is still near to the value 2 separating region 2 from region 1 (cf. Fig.|3J. 
Second, for values of K/D around the transition from region 2 to region 3 the system apparently relaxes very slowly 
to its asymptotic behaviour and it becomes very difficult to extract reliable values for (Nz) from the numerics. The 
dependence of (Nz) on the ratio K/D as obtained numerically is summarized in Fig. El 



6 




FIG. 2: Typical time dependence of the collective orientation S in the three different regimes described in the text. Shown 
are Sx(t) (dashed lines) and Sy{t) (solid lines). The parameter values are = 0, a = /3 = and K = 3. (a) Re eime 1: 

D = 2.4, (N^) = (b) Regime 2: D = 1.2, (N^) = 3.5 10"^ (c) Regime 3: D = 0.6 {N^} = 0. 




FIG. 3: Numerical results for the time averaged torque {N^} for = as a function of the disorder strength D tor K — 3,a — 
(3 = l/y/2,H^ = 0. The vertical dashed lines divide the different regimes from each other (1.5 < D, regime 1; 0.8 < D < 1.5, 
regime 2; D < 0.8, regime 3). In the interval indicated by the dotted lines, 0.75 < D < 0.86, no asymptotic value of (A^'^) could 
be obtained after 200 periods of driving. 
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IV. EFFECTIVE FIELD APPROXIMATION 



In order to obtain an improved understanding of the numerical results obtained in the previous section we now turn 
to an approximate analytical analysis of the system. A valuable tool in this respect is the effective field approximation 
for ferrofiuids as extensively reviewed in To apply this approximation in the present context the following steps 
are performed. From the central equation Ijl9ll we derive an exact equation for the time evolution of (u) of the form 



dt {u) +D{u)=- {u{u ■ H,nf)) + H, 



mf ■ 



(34) 



where 'H.-mf is determined by l|2(11l . In order to decouple the higher moment on the r.h.s. of l(3ijl we approximate the 
unknown distribution P{(f>, t) by an instantaneous equilibrium distribution of the form (|21ll corresponding to a so far 
undetermined effective magnetic field He(t) 



P(0,t) 



1 



2vr/o(^) 



cxp 



ge,x cos(0) + He,y sin(0) 
D 



(35) 



Using this distribution to calculate the averages in l(3ijl and using (cf. l|22 



D J 

the following evolution equation for the effective field He(t) may be derived 

H, 



dt 



B ^ 



He 



^ ^^H, X (He X U^f) ^ ^ 



(He — Hm/) 



With the help of 



H 



mf 



D 



(36) 



(37) 



(38) 



which follows from ll271|l and ll36|l and representing He by modulus and phase according to He = He{coBip,sa\ip) we 
find the following closed set of nonlinear ordinary differential equations for the time evolution of the effective field 



D 



' He 



— B — 



D 



D 



-B 



El 
D 



B' 



He 

D 



- KB 



D 



Hy (t) sin (p 



D 



Hy{t) COS (fi 



(39) 



(40) 



These equations cannot be solved analytically, however their numerical solution is much simpler than the numerical 
solution of the Fokker-Planck equation Ijlijll . As shown in Fig. 01 the results provide rather accurate estimates for the 
relevant quantities. 



V. DIFFERENT MANIFESTATIONS OF THE RATCHET EFFECT 



With the help of the effective field approximation it is possible to gain a more intuitive understanding of the 
existence of different regimes of the system behaviour as found in section ITTTI Moreover, their stability as well as the 
operation of the ratchet effect in the different regimes may be elucidated. 

Using the evolution equation for the effective field II39|I . H4()II and the relation Ij36|l between the average collective 
orientation S and the effective field we may write the time evolution of S in the form 

dtS = -yViS)+Fili,S) . (41) 
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0.5 



1.5 



FIG. 4: Time averaged values of the a;-(solid lines and crosses) and y-component (dashed lines and circles) of the average 
orientation S as function of the noise strength D. The other parameter values are as in Fig. [2| Lines are numerically exact 
results from the Fokker-Planck equation 11!) |l whereas symbols denote the results obtained from the numerical solution of the 
effective field equations l|39|l and l|4f)|l. 



The r.h.s. of H41II was split into a central potential term incorporating the effects of diffusion and interaction and an 
external force field related to H(t). Note that < S' < 1 always. 
Let us first discuss the potential part. It is given by 

n^)-f^^-A'f (42) 



where B^^ denotes the inverse function of B defined in l|23ll . Hence B~^ is linear for small values of S and tends to 
infinity for 5 — * 1 (cf. lf2ijl . lf25ll ) . For small values of K/D the potential V{S) therefore has a minimum at S* = 
whereas for K/ D > 2 it develops a non-trivial minimum at Seq > (cf. Fig. left column) . From Ij38|l and Ij4()|l the 
corresponding equilibrium value of the effective field satisfying 



is determined by 



B[^]^S,q (43) 



H^^^KBl^]. (44) 



D 



We note that with increasing K also Seq increases (see Fig. El . In the absence of an external magnetic field, H — 0, 
the stationary solution Seq of the effective field equations coincides with the exact equihbrium solution of the Fokker- 
Planck equation Ijlijll . 

For the external force field we get 

F{t) - -^4^Hy(*) e, + (^1 - -^^j^^ Hy{t) cosy, . (45) 

For small values of S this field is parallel to the external magnetic field H. For 5* 1 on the other hand the radial 
component in F is suppressed and only the tangential term survives (cf. Fig. El middle column). This is quite intuitive 
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X X 



FIG. 5: Plot of the different terms contributing to the r.h.s. of l|41|l for Hx = K = 2> and a time t for which Hy(t) — 1. The 
first line corresponds to D = 1.6 (regime 1), the second one to D = 1.2 (regime 2), and the third one to D = 0.6 (regime 3). 
Shown are the potential 1^(5*) defined in Il42|l in the left column, the external force field F defined in Ij42|l in the middle column 
and the total resulting force field in the right column. The dashed line in the middle column marks the equilibrium value Seq- 

since for S' ^ 1 all magnetic moments are aligned and therefore the external field induces identical changes to their 
orientations giving rise to a change in the orientation of S only. 

For small values of the external magnetic field S will only slightly differ from its equilibrium value Seq. We will 
hence find 5 ~ for small values of K/D which was characteristic of regime 1. Larger values of K/D result in values 
of S differing substantially from zero. From it is then clear, that = 7r/2 is a possible steady state since in 
this case F is parallel to and there is hence no systematic force defiecting S from the orientation along the y-axis. 
This situation corresponds to regime 3. Finally, there is also a steady state solution with s > and (/s ~ (regime 2), 
which is, however, not obvious from I^SJ- 

We will now give a simple argument for the transition between regime 2 and regime 3 which allows to derive an 
estimate for the threshold value of K/D separating these two regimes. Since the potential part in II41II only infiuences 
the modulus of S the time evolution of the angle (p is solely given by 

dtf = . (46) 

Assume now that at some given point of time we have < ip{t) < 7r/2. As long as Hy{t) is positive (as in Fig.l^l S 
will be pushed in positive y-direction (cf. Fig.0 right column) implying a slight increase of S beyond the value Seq- 
Similarly, when somewhat later Hy(t) changes sign and becomes negative S will be pushed in negative y-direction 
and S will slightly decrease. Simultaneously ip changes according to (j46|l . i.e. increases when Hy{t) > and decreases 
when Hy{t) < 0. However, the rate of change depends on S{t). A rough estimate on whether the increase of ip 
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dominates over its decrease or vice versa may be obtained from the quantity 



Jiip) 



■ V 



dS \Hy{t)S 



(47) 



If J{ip) < the decrease of if will dominate, if J{ifi) > its increase. Remarkably, although J depends on ip its sign 
depends only on the ratio K/D. Hence either tp decreases all the time until it reaches (p = and we are in regime 2 
or it increases up to ip = tt/2 corresponding to regime 3. An explicit calculation yields 

= H!^{K^ - DK - H!^) ^"^^^ 

As shown in the appendix, the denominator is positive for all values K/D > 2. The condition for the stability of 
regime 2 J < is hence 

+ KD + _ < , (49) 
which together with H44II yields that regime 2 is stable if 



, K \ D 



or equivalently 



^<3.75.... (51) 



As we will see in subsections IV Bl and IV ("I below a more systematic investigation of the stability of regimes 2 and 
3 will give rise to the same result (cf. l(85ll . lllfl0ll V 



A. Regime 1 

Regime 1, characterized by K/D < 2, is the most relevant for ferrofluids. There is no spontaneous collective 
orientation, i.e. S = 0, and accordingly we must have Hx 7^ in order to find a noise induced rotation. The situation 
is hence similar to the case without interaction, nevertheless we will show that the interactions bring about a strongly 
reinforced ratchet effect. 

To obtain a quantitative estimate of this reinforcement we use a variant of the perturbation theory for small values 
of the external field introduced in . To keep track of the different orders in the expansion it is convenient to use 
the rescaling 

H ^ eH . (52) 
We then solve the Fokker-Planck equation (tT9|l using the expansion 

P = P("' + eP(i) + + £3p(3) + ... (53) 

with the unperturbed solution given by 

Using l(27|l this ansatz gives rise to a similar expansion for the coefficients a„ with ?i > 

an = ealP + e^a^f^ + ... (55) 
whereas bq is to all orders in e fixed by the normalization condition to 

ao^^. (56) 
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The pecuHarity of the first regime is that S ^ when H 0. We may therefore consistently employ a similar 
expansion for the collective orientation S, 

S ^ eS^ +e^S^ + ... . (57) 



The first non-zero result for (A^^) is obtained in fourth order in e. From we infer that we hence need ai up to 
third order. Similar to Q we find to first order in e 

a + Z,-f)„».i^ipi, ,58, 



and to second order 

(^dt + D- y j aY> = (59) 

{dt + 4:D) ^ {H^ ~ iHy)a[^^ + 2TTKa[^^a[^^ (60) 
In the first regime we have D — K/2 > 0. From Ij59|l we then get a''^\t) for t oo. To third order we find 

a , n ^'^\ (3) H^+iHy (2) (1) (2) , . 

dt + D^ — ]a\'= ^4 - 7rXa!_i4 . (61) 



2 ^ 

Using the special time dependence ^ and solving for the asymptotic behaviour of a^^^"* and ap'' we find for the torque 
(IHOll after some algebra 

{128D^K - QAD'^K'^ + iOD^K - 20D^K^ + SDK - 2K^) sin S 

(62) 



(161)2 + 1)(4D2 _ 4£)^' _^ ^ 4)2(2£) _ K){4:D^ - ADK + + 16)(4D2 ^ 
184i:>5 ^ QgD^K + 58D^K^ + 200D'^ + 8D^K - 3D^K'^ + 16D + 6DK^ + iK) cos (5 1 



(161)2 _^ i)(4£)2 _ 4£)x + X2 + 4)2(2D - X)(4Z)2 - ADK + K^ + 16)(4D2 ^ 
For i<r = this expression simplifies to j3| 

v2 fl (O'i n2 



3F:,a2^ (23L»2^_2)cos(5 



^^'^ 8 {4D^ + 1){D^ + 4){D^ + 1){16D^ + 1) ' ^^^^ 

which coincides with the result found in when specialized to the x-y plane, i.e., to 9 = -k 12. Fig.Elshows the torque 
as function of the phase angle b for different values of the interaction strength K. It is clearly seen that although 
the interaction is in regime 1 too weak to induce qualitative changes in the behaviour the value of the noise-induced 
torque is greatly enhanced by the interaction between the particles. 



B. Regime 2 

In regime 2 the interactions between the particles give rise to a non-zero value of the x-component of the 
collective orientation S. This breaks the symmetry x i-> —x even in the absence of an x-component of the external 
field, = 0, and results in an interaction sustained ratchet effect in the system. In the present subsection we derive 
an approximate solution of the effective field equations (I39|) . lj4()|l characteristic for this regime and discuss its stability. 

To obtain an approximate analytical solution valid in the long-time limit we resort again to perturbation theory and 
assume that H is much smaller than the equilibrium value H^q of the effective field He- This assumption can easily 
be met if we are sufficiently far from the boundary of regime 2, K/D = 2, at which Heq tends to zero. To formally 
organize the perturbation expansion it is again convenient to employ the rescaling (I52|l and to write the solution of 
the effective field equations Ij39|) . Ij40ll as power series in e 



i7e = i/i°)+ei?i^)+e'i/P + ... 



(64) 
(65) 
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FIG. 6: The time averaged torque in regime 1 as function of the phase S in the time dependence ^ for different values of the 
coupling strength K. The symbols are results from the numerical solution of the Fokker-Planck equation the lines show 

the results of the perturbation theory outlined above. The parameter values are ^ a = P = l/Ve, D = 1.6, and /f = 
(dashed line and circles), K = 0.5 (solid line and crosses) and K = 1 (dashed-dotted line and diamonds) respectively. 



where = Heq and (/s^"^ is not specified at this stage. Using these ansatze in . l|171|i and matching powers of e 



eq 

we find to first order 



n\ — — 2KD „ sinwn 



a*™ = '^"'^'"-'"'V .w ^ (67) 

J^eq 

These equations can be solved for arbitrary (po ■ The second order equations are rather long and will not be displayed 
in full generality. However, as typical for degenerate perturbation theory, the second order equations may contain 
secular terms. In the present case this happens unless either ipo = ±tt/2 or J^^ dt Hy{t)He^\t) = 0. The first 
condition corresponds to regime 3 to be discussed in the next subsection. The second one is equivalent to (/jq = as 
follows from and is therefore the condition appropriate for the present investigation of regime 2. Using i^o = 
and the specific form |(2J of the time dependence of Hy(t) the first order solution tends for large t to 

i?^^^ = (68) 
(1) K-D f /? 



V^'' = [a sin(t) - I cos(2t) J + , (69) 

where is an integration constant. Using these results the second order equations greatly simplify and take the 
form 



dtp^^^ = , 



(70) 
(71) 
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with 

K"^ - 2KD - Hi, 
^--Hl^ + KD-K^- 

Hence ip^'^'^ is a constant and i/i^^ may easily be calculated for all values of Lpm- However, this solution produces 
secular terms in the third order equations unless 

r dtHy{t)Hi'\t)^Q (74) 
Jo 

implying 

^ 3/3a^(A^ + 2) K^D 

This result for (pm completes the first order solution of the effective field equations. Higher orders may be obtained 
along the same lines but will not be considered here. 

For the collective orientation we then obtain to first order in the external field 

S. = ^ (76) 

Sy - ^^^'^ ■ (77) 

From we hence infer that the time average of Sy is directly related to (pm , 



^ El 
K 



Sy = ^<Pm. (78) 



Fig. 13 compares the first order results obtained above with numerically exact results from the Fokker-Planck equation 
ifToll for different parameter values. The agreement is rather good showing that both parameter sets for a and (i 
correspond to external field strengths well within the region of validity of the perturbation theory. 

We now turn to the investigation of the stability of the first order solution given by Ij68|l . lj69|l . and Ij75|l . To this end 
we add small deviations 5He and Sip to this solution and study their time evolution by Hnearizing (l39)l . (liHI around 
Hi^\t),ip^^\t) [Hi. In this way we find 

(5H,\ _ ( an -auHy{t)\ f SHA 



where 



H'l, + 2DK - 

eq 



KD + Hi, + - /\ 2 
a.. - j}^ . (82) 

We are interested in the stability of the fixed point SHe = S(p = oi this system. Because of the periodic time 
dependence of the coefficient matrix in iTTQ)! a general analysis requires Floquet theory ^|. However, the special form 
of H79|l allows us to use a more direct method. 

We first note that one can show (see appendix) an < within regime 2. We next define the function 

V{SHe,5p) = ^{a2iSH^ + aiaV) ■ (83) 
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FIG. 7: (a) Time average of the y-component Sy of the collective orientation S in regime 2 as function of the noise strength 
D for interaction parameter K — 3. Symbols are numerical results from the Fokker-Planck equation, lines show perturbative 
results from the effective field equations. The solid line and the crosses are for a — (3 = 1/V2, the dashed line and the circles 
are for a = P — 0.3. (b) Time evolution of the y-component of the orientation vector S in regime 2. The noise intensity is 
D = 1.2, all other parameter values are as in (a). 



For its time derivative we find 



Moreover if 



we have for all values of SH^ and Sip 



^ = al.anH^ < . (84) 



021012 > (85) 



V{6He,Sip)>0 (86) 

with equality being valid only for SHe ~ 5ip = 0. V{SHe,5ip) is hence a Lyapunov function j23| of the system ifTfljl 
and the point SHe = Sip = is stable as long as l(85|l holds. From (IsTll one finds ai2 < (see appendix) within regime 
2 and hence the first order solution He^\t) , ip^^'^ (t) is stable if 021 < 0. From II82|I it is seen that this condition is 
identical with Ij49ll . We hence find again that regime 2 is stable if 

2 < ^ < 3.75 . (87) 
Note that the entire stability analysis performed above is valid only to first order in the external field strength. 



C. Regime 3 



Similarly to the previous subsection we are now using the effective field equations Ij^ffl|l . lj4()|l to derive an approximate 
analytical solution typical for regime 3 and to discuss its stability. In regime 3 we have ip = n/2 (the case ip = —Tr/2 
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0.92 




FIG. 8: Time evolution of the y component of the orientation vector S in regime 3. Symbols are numerical results from the 
Fokker-Planck equation, lines show perturbative results from the effective field equations. The parameters are D — 0.6 and 
-ft" = 3, as well as a = /3 = 0.3 (crosses and full line) and a = f5 = 0.1 (circles and dashed line)respectively. 



can be dealt with similarly) and we find from Ij39|l . (14011 



' He 



dtip = 



\ D 



\ D \ D 



H, - KB 



He 
D 



Hy{t) 



(89) 



From the second equation we get immediately Lp = 7r/2. The first equation determines He{t) and cannot be solved in 
closed form. Focusing again on the situation of small external fields we use (I(i4|l and obtain for the first order term 
the equation 



dtHi^'^ = -AHi^^ +BHy{t) 
with A and B defined by (I72|) and Ij73|l respectively. The solution is given by 

Hi^^ =Be-^' j dtHy{t) e^* . 

For the special time dependence Q we find for large t 

^(1) ^ ^^ Acos(t)+sin(t) ^ ^^Asin(2t) - cos(20 



+ 1 



^2 



(90) 



(91) 



(92) 



In Fig. ISI this result is compared with the solution of the Fokker-Planck equation showing again good agreement for 
the parameter values chosen. 

We now turn to the investigation of the stability of the first order solution. In view of the form of the potential 
V{S) in regime 3 (see Fig. |21 first column, last line) it is reasonable to assume that the state will be destabilized by 
perturbations in ip. Restricting ourselves to this case we use He{t) as determined above, set ip{t) = 7r/2 + Sip{t), and 
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determine the linearized time evolution of S(p{t). From (j4()|l we get 



with solution 



D fHe 
1 B — 

dtSv = Hy{t) 6^ (93) 



D n (H, 



Sifit) = 5^0 exp dt' ^p^Hyj . (94) 

Since both Hy and are 27r-periodic functions of time this solution may be written in the form 

S^{t) = S^o e--^"* fit) , (95) 

where /(i) is also 27r-periodic and Cq is the constant term in the Fourier expansion of the exponent in and 
therefore given by 

ao = -y^ dt ^^^^^ Hy. (96) 
Using the first order solution He,{t) ~ Heq + H'^\t) we find 

'''\t' Hi^\t') Hy{t') . (97) 



ao 



The integral in this expression is in regime 2 and in regime 3 always positive as can be shown by expanding Hy{t) in 
a Fourier series 

oo 

and using l(9T|) . We then find 



n— — oo 



^ |2 



1 r ' — ' 

-/ dtHi'Ht)Hyit)^j2 

^0 -V,— n 



n— 



|6„r > > (99) 



since B > as demonstrated in the appendix. Hence the sign of oq as given by l(97jl depends solely on the prefactor 
and the solution is unstable if 

DK + - + Hl^ < . (100) 

This condition is again identical with if^Cijl and we hence find that regime 3 becomes unstable at the same value of 
K/D at which regime 2 becomes stable. 



VI. CONCLUSIONS 



In the present paper we have theoretically analyzed the infiuence of particle-particle interactions of the mean-field 
type on the ratchet effect in ferrofiuids. We have used the simplified mean-field model because the details of the 
realistic dipole-dipole and hydrodynamic interactions are too complicated for a general discussion. On the other 
hand, several qualitative effects which may be expected to show up in case of these realistic interactions can already 
be discussed within the mean-field approach. 

Quite generally we have found that interactions that favour the parallel alignment of the magnetizations of the 
ferrofiuid particles reinforce the ratchet effect. In the present situation this means that the magnetic torque per 
volume of the ferrofiuid is enhanced. Although intuitive this result in by no means trivial since the coupHng between 
the ferromagnetic grains is also likely to reduce their orientational fiuctuations which are the driving force of the 
ratchet effect. 
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Moreover we have shown that for sufficiently strong coupUng giving rise to a spontaneous symmetry breaking with 
respect to the collective orientation of the grains the ratchet effect may be brought about for situations in which it 
could not operate in the absence of interactions. In the present system this happens for a purely oscillating external 
magnetic field without constant component in the a;-direction. In this case and without interactions no ratchet effect is 
possible due to symmetry reasons whereas a spontaneous breaking of the relevant symmetry caused by the interactions 
may induce a rectification of fiuctuations. 

For even stronger coupHng strength between the particles this self-sustained ratchet effect again disappears. The 
particle orientations are now highly aligned and closely follow the direction of the external field. Fluctuations are 
therefore suppressed and no rectification is possible anymore. 

Our analysis builds on the non-linear Fokker-Planck equation lfT9|l for the collective orientation (tT5|l of the ferrofiuid 
particles. A first insight into the behaviour of the system is gained from a numerical solution of this equation. We then 
use a variant of the effective field method, which is a well-known tool in the theory of ferrofiuids, to get approximate 
analytical results for the transitions between and stability of the different regimes of operation of the ferrofiuid ratchet. 

The central parameter distinguishing the different regimes is the ratio between the dimensionless coupHng constant 
K and the intensity of the fiuctuations D. For K/D < 2 no self-sustained ratchet effect is possible since the interactions 
are too weak to induce a spontaneous symmetry breaking. For 2 < K/D < 3.75 a self-sustained ratchet effect may be 
observed with a wide spectrum of values for the transferred angular momentum. For 3.75 < K/D the self-sustained 
ratchet effect disappears again due to an instability in the steady state solution for the collective orientation of the 
particles. 
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APPENDIX A 



In this appendix we determine the signs of some expressions needed to determine the stability of regimes 2 and 3. 
In these regimes we have K/D > 2 the value of Seq is determined by the non-trivial minimum of the potential V{S). 
From (I42|l we find by differentiation 

verifying that Seq depends only on the ratio K/D. Using the properties of the function B{x) defined in lf23jl one can 
show that for all S G [0,1] 

This in turn implies 

K^ - DK - K^Seq > (A3) 

and using l(i3jl we find 

K^ - DK - Hl^ > . (A4) 

Due to Ij73p and Il8ip this immediately implies 

^--^-^-- H^^ + KD^K^ ^'- 
Using again the properties of the function B{x) one can also show 

l-2^<^^ (A6) 

Similar manipulations as used above then yield for K/D >2 

K^ - 2DK - Hl^ < (A7) 
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from which by using H72II and H8()|l we find 



DiK"^ - 2KD - Hi,) 
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